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Figure 1: Diffractive vector meson photoproduction at large momentum transfer. 
1. Introduction 

In a previous paper, |||], we performed a detailed theoretical study of the diffractive produc- 
tion of vector mesons at high momentum transfer in 7p collisions. We worked throughout 
in the leading logarithmic BFKL framework and factorised the meson production from 
the hard subprocess using a set of meson light-cone wavefunctions. 

The process is illustrated in Figure The photon and proton collide to produce 
a final state containing a vector meson and the products of a proton dissociation. The 
meson is produced with large transverse momentum, the square of which is equal to the 
Mandelstam variable —t, and this is in turn much smaller than the available centre-of- 
mass energy s, i.e. s ^ —t ^ Aq^j-,. This hierarchy ensures that the meson and the proton 
dissociation are far apart in rapidity. Note, since we will always integrate over the products 
of the proton dissociation we do not exclude the possibility that the proton may scatter 
elastically. However, this contribution is negligibly small for sufficiently large t whence 
the finally state is dominated by configurations containing a single jet with transverse 
momentum balancing that of the meson [^. 

This process has been subject to considerable experimental investigation in recent 
years, with the measurement of the t-distribution of the meson (p, (j) and J/^) |, |] and 
the spin-density matrix elements extracted from its decay |Q, ^ being the highlights. It is 
the purpose of this paper to compare the theoretical results of |||] to the data. Before doing 
so, we recall the status of diffractive meson production at high-t prior to this analysis. 

The experimental data indicate that the cross-section da/dt ~ for p and cj) produc- 
tion. Moreover, the relative smallness of the Tqq spin-density matrix element suggests that 
the interactions tend to produce transversely polarised mesons. This is in stark contrast 
with simple theoretical expectations. Lowest order perturbative QCD does indeed predict 
a ~ distribution but for longitudinally polarised mesons. For transversely polarised 
mesons, the expected distribution is ~ This observation has led the authors of Q to 
postulate that the production of transverse mesons is enhanced by the presence of a large 
non-perturbative coupling of the photon to chiral odd quark-antiquark configurations^. 

^We refer to "chiral even" and "chiral odd" configurations to denote the coupling of the photon to a 
quark and antiquark of equal (chiral even) and opposite (chiral odd) chirality. The chiral odd coupling 
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The authors of |^ also identified the breakdown of perturbative factorisation in the two- 
gluon exchange model due to divergent contributions from configurations where either the 
quark or antiquark carry all the momentum of the meson (i.e. end-point contributions). 
They argued that this divergent behaviour would be tamed by Sudakov effects and that the 
dominant chiral odd amplitude was in any case free of problems. In contrast, the authors 
of make a qualitative argument that it is precisely these end-point contributions which 
are responsible for the behaviour seen in data rather than the which would be 
expected in perturbation theory. Finally, we should also comment on the apparently naive 
approach of In this analysis, longitudinal meson production is suppressed by virtue 
of the fact that the quark and antiquark share equally the longitudinal momentum of the 
meson, and transverse meson production is enhanced due to the use of the constituent 
quark mass. Consequently, in what follows we shall take care to comment on each of the 
above analyses as appropriate and we shall pay particularly close attention to the end-point 
behaviour; there being no end-point divergences in the BFKL treatment, even for massless 
quarks. 

Before proceeding we should summarise the experimental analyses. The ZEUS data 
are for the energy range 80 GeV< W-^p < 100 GeV with the mean value approximately 
central. Their measurements are for 1 GeV^ < — t < 10 GeV'^ and they make the cut 
-(M|. - ml)/t < 100 where Mx is the mass of the proton dissociation products and rup 
is the mass of the proton. The relation 

Ml-ml = -t{l/x - 1) (1.1) 

translates this to a cut of x > 0.01 where x is the longitudinal momentum fraction of the 
struck parton. ZEUS quotes a photon virtuahty Q"^ < 0.02 GeV^ with (Q^) ^ iq~5 GeV^. 
This allows us to neglect in the kinematic relations and in our theoretical calculation. 

The HI data are for J/ip production, and in the energy range 50 GeV < W^p < 150 GeV 
with the mean value approximately central. Their measurements are for 2 GeV^ < —t < 
30GeV^ and they cut in the 'elasticity' variable Ze, which is related to Mx to a good 
approximation (for the J/ip) by 

Ze^l-(Mi-t)/1^4. (1.2) 

HI quote Mx < 30 GeV which translates to Zg > 0.95 and x > —t/{l — Zf.)s. The photon 
virtuality is restricted to satisfy <1 GeV^ with (Q^) ~ 0.06 GeV^, again small enough 
to regard the photon as being real. 

Following , we write the cross-section as a product of the parton level cross-section 
and parton distribution functions: 

da{-ip VX) ( r^f ^ ( - , ^^\ ^^(T? ^ ^<l) 



dt dx 



[j^^^G{x,t)+Y,Mx,t)+qf{x,t)]^ ^^^^y ^'^^ ' (1-3) 



where Nc = 3, G{x,t) and qf{x,t) are the gluon and quark distribution functions respec- 
tively and s is the jp centre-of-mass energy squared. The struck parton in the proton, that 

vanishes for a pointlike coupling to massless quarks. 
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initiates a jet in the proton hemisphere, carries a fraction x of the longitudinal momentum 
of the incoming proton. 

The parton level cross-section, characterised by the invariant collision energy squared 
s = xs, is expressed in terms of the helicity amplitudes Mw^s, t) (A is the photon helicity 
and A' is the meson helicity): 



1 



da 



dt 167rs2 



(|M++(s,t)|2 + |M+o(s,t)|' + |M+_(s,t)|2) 



(1.4) 



The experimenters also measure the spin density matrix elements r^j' which are extracted 
from the angular distribution of the decay products of the vector meson: 



d^a 



oc 



i(lTrOo^)±i(3rOo^-l)cos2 0, 



d cos 9h d(j)h 

=F\/2Re[roo] sm29h cos =F r^i^ sin^ 9h cos2(j)h 



(1.5) 



where the upper (lower) signs are for spin-0 (spin-1/2) particles. 6^ and are the spherical 
polar angles of the positive particle of the two body decay of the meson, where the vector 
meson momentum defines the z-axis and the x-axis (which fixes (ph = 0) is defined to lie 
in the direction of the hard momentum transfer, q. Note that |qp = —t in what follows. 
The r-matrix elements can be written 
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(1.6) 
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(1.8) 



where (...) denotes the integration of the parton level quantities over partonic x with the 
appropriate cuts. Throughout this paper we choose a fixed value for the 'yp centre-of-mass 
energy: = W^p = 100 GeV. 



2. Theoretical Results 

We begin by summarising the results of [|l| for the relevant matrix elements. We append a 
superscript to specifiy if the photon-quark coupling is chiral even or chiral odd. 

We use the meson distribution amplitudes of [12,|lS] and use their notation throughout. 



The explicit forms for the distribution amplitudes we use can be found in Appendix A. In 
what follows we systematically include terms of progressively higher twist in the distribution 
amplitudes in order to ascertain their relative importance. The relevant amplitudes are 
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n=+oo oo u'^ + n^ exp[x2n(t^)^] ^ / o n 

1^^ ^^[z.2 + („_i/2)2][^2 + („ + i/2)2] ^0-1 2n, 1), 

(2.1) 

= '-^^ [ du - (1 - 2u) [ Mh (-) - (-))) 
"V^~ r , ^^ + »^ exp[x2n(^^)^] ... 



n=— cx) 



(2.2) 



^ ^ / '^"[.2 + (n-l/2)2][.2 + (n + l/2)2] sin(z..) ^n, 1). 

(2.3) 

d^;(/i[*^(i;)-<^x(^)) 

+ e> 
2 + (n- 1/2)2] + (n + 1/2)2] ~sin(i7rz/) -f-f 



n^oo ^2 ^^2 exp[x2n(^)^] , 



(2.4) 



i^^ + n2 exp[x2>l^] 

'2 + (n - 1/2)2] [z/2 + („ + 1/2)2] sin(i7rz/) 



(2.5) 



n=+oo »oo ,,2 I ^2 



X 



r , '^' + n' exp[x2n(^-)^] r f o m 

£"^i-oo K + 1/2)2] [z.2 + (n + 1/2)2] ^^^^^^^^ l3_i[i^,Zn,q,u,U). 

(2.6) 
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Qv 


1/V2 
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2/3 



Table 1: The meson electric charge couplings. 



where 

I^p{u,n,q,u;a) = mj <f p p^+^ p* Ka{m\p\)e'i^''' P-^'^"*^ 

X [J^(gV/4)J,,(<7p74) - (-1)" J_^(gV/4)J_^(gp74)], (2.7) 

= 2u — 1, p = n/2 — iu, p = —n/2 — iu, Ka{x) is the modified Bessel function and the 
parameter a equals 1 for the chiral-even and for the chiral-odd contributions. In ( |2.7D we 
use complex variable notation for the transverse vectors p and q, e.g. p = Px + ipy 
After some effort Q one finds 

n, q,u;a) = - ^r{a/2 - C)r(-a/2 - C) {i sign (1 - 2n))°-/^+- 

X ^1^^ [sin7r(a + ^ + C) B{a, p,q* ,u,C) B{P,p,q,u* ,C) 

- {-irsimria- p + C) B{a,-p,q*,uX)B{(3,-p,q,u*,C)] (2.8) 

where we have introduced the dimensionless parameter Tq = ^m} l\q^ and the conformal 
blocks 

^(a,,,,*,.,C) = (-4n.)-(^— 0/. (^)°2-- "^^ + ^ + y^^ 

^^^i 2 ' 2 '^ + ^'4^J- ^'-'^ 

Note, that the sums are performed over even conformal spins 2n. 
We have also introduced the energy variable 

3a= , / 2;s \ , , 

where A is an undetermined energy scale in the leading logarithmic approximation. The 
eigenvalues of the BFKL kernel are proportional to 

Xn{y) = 47^e(^V(l)-V'(l/2 + |n|/2 + iz.)^ (2.11) 

and we collect together several factors into Cy. 

iV^ - 1 

Cv = ial^^eQv (2.12) 

where Qy is determined by the electric charge of the quarks which make up the meson. 
The explicit values are listed in Table ||. 



- 5 - 





P 


<A 




fv [GeV] 


0.216 


0.231 


0.405 


/^(IGeV) [GeV] 


0.160 


0.215 


(0.405) 



Table 2: The meson decay constants ]l2[ . 
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no- flip 


j^odd 


j^even 




double-flip 







^odd 



Table 3: Classification of the relevant hclicity amplitudes in twist. 



All the helicity amplitudes depend upon the vector meson coupling constants; either 
fv (vector) or fy (tensor). Table ^ presents the vector and tensor couplings for the relevant 
mesons. The vector couplings are directly measurable from the electronic decay width of 
the vector mesons whilst the tensor couplings have been calculated at ^ = 1 GeV using 



QCD sum rules |ll2| 



2 



In Table ^ we classify the above set of helicity amplitudes in terms of twist. We work 
to next-to-leading twist in each amplitude, i.e. we do not consider twist-4 terms in the 
single and no flip amplitudes. 

3. Parameters and prescriptions 

Strictly, in perturbation theory, we should use the current quark mass. However, using the 
current mass renders the chiral odd contribution negligible for light mesons and proves to 
be incapable of describing the r-matrix elements. The authors of noted the importance 
of introducing a large chiral odd contribution. The method they employ introduces a large 
non-perturbative coupling of the photon to chiral odd quark-antiquark configurations. We 
shall take a different approach and replace the current quark mass with the constituent 
quark mass, i.e. we take m = my/2. We will investigate the sensitivity to this parameter. 

We also need to explain how we implement the strong coupling, and the scale A'^ 
which appears in ( p.lOf ). We might reasonably expect that A^ = Pmy — where (3 and 7 
are unknown. To simplify matters, we choose to fix, /? = 1 since a change in its value can 
approximately be traded off with a change in q^. 

Notice that Og appears in two places. It appears as a prefactor of in every helicity 
amplitude. It also appears in the definition of z (see (^l0|) ). The prefactor arises from the 



coupling of the two gluons to each impact factor. The coupling in z, however, is generated 
by the gluon couplings inside the gluon ladder. We shall denote these couplings al^ and 
^BFKL respectively. Strictly speaking a fixed value of as is appropriate to the level of our 



^The tensor coupling for the J/ip value is unavailable in the literature and so we, as default, put it equal 
to the vector coupling. 
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calculation. However, next-to-leading logarithms cause ag to run. Arguments exist that 
the effect of these corrections is to effectively fix af^^^ [14|. Phenomenologically, fixing 



^BFKL j^g^g pj-Qved successful in the past 15| and we will fix it here. We remain free 
to run or fix the impact factor coupling, al^ , and we shall subsequently make use of this 
possibility. When we refer to a running coupling in future, we shall be referring to its value 
at 1 GeV. We note that, as ratios, the r-matrix elements are independent of al^ . 
In summary, the parameters we shall treat as free for fitting purposes are 

ai^, af^^^ and = - 7*. (3.1) 

We shall also show the sensitivity of our variables to the quark mass m though we shall 
not treat this as a fit parameter. 

For other parameters, we shall use either those measured by experiment or estimated 
in the literature. We shall not vary these. The sensitivity of our predictions to various 
approximations for the meson distribution amplitudes will be tested by testing the impor- 
tance of l/|qp suppressed and above terms and by employing three different distribution 
amplitude prescriptions. We refer to the three prescriptions as the 5-function prescription, 
the asymptotic distribution amplitude prescription (obtained in the limit — > 0) and 
the full distribution amplitude prescription (of |l3|). All necessary information on the 
meson distribution amplitudes can be found in Appendix |^. 



4. Phenomenology for the p meson 

In this section we focus wholly on the phenomenology of the p meson. We begin with the 
most crude approximations for the meson wavefunction and systematically relax these to 
the point where we eventually treat the meson using the full light-cone wavefunctions of 
H, PI- 



4.1 The collinear and 5-function approximations 

The collinear approximation corresponds to putting the transverse momentum of the quark 
and antiquark exiting the hard scatter (measured relative to the meson direction) to zero. 
This approximation generates J-functions which force the transverse size of the quark- 
antiquark pair to zero. This fact forces all the higher twist helicity amplitudes to zero so 
that only M^^''" and M""^ (twist-2) then survive. 

The (5-function approximation for the meson distribution amplitudes corresponds to 
equating 

(t>\\{u) = Mu) = Hu-l/2). (4.1) 
The MV'q^ helicity amplitude is zero in this approximation and so the only non-zero am- 
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plitude is M°^: 



Ml^^ = sCv^ fvmv 



47 J {u'^ + {m-l/2Y){v'^ + {m + l/2f] 



ds ^i/2+s+\m\ r(l - s- iv) T{1- s + iu) 



■ r 



,^2m T{1- s/2-iu/2)T{l- s/2 + iv/2) 

r^(l/2 + s + |m|) 

r(l/2 + s/2 - iv/2 + |m|) r(l/2 + s/2 + iv/2 + |m|) ' 



(4.2) 



Figure |2| shows the predictions for da/dt for three different sets of parameter values. 
Running the scale that appears under s in the BFKL logarithm ( 2.1C| ) with t or running 



al^ with t prove to predict t dependencies that are too steep and incompatible with the 
data presented by ZEUS [Q]. Fixing the strong coupling in the impact factor and in the 
BFKL ladder in addition to the BFKL logarithm yields a good fit. These results agree 
with those of Q for the leading conformal spin component (n = 0)^ and [^] for n > 0. 
Since only the M?*^ amplitude is non-zero with the collinear and (5-function approxi- 



mations imposed, referring to their definitions (( |1.6[ ), (1.7), (|1.8D ), we see that this enforces 



all the r-matrix elements to zero. Figure § shows this prediction in comparison to the 

00' 



ZEUS data. We see that zero r-matrix elements are compatible with data for Tqq, in fair 



agreement for Tj^q and in poor agreement for r^^^i- 

4.2 Collinear approximation and asymptotic distribution amplitudes 

The use of a J-function distribution amplitude is a crude approximation. For the charm 
quark and J/ip meson, rric ~ mj/^/2. For this meson, the configurations that provide 
the dominant contributions should be close to n ~ 1/2, i.e. we expect the (5-function 
approximation to be fair. This is not so for the light p. The use of the (5-function was 
motivated purely by past phenomenological success To understand this success, and to 
have any hope of predicting the helicity structure of the cross-section, we must relax this 
crude assumption. For now, we maintain the collinear assumption and use the asymptotic 
distribution amplitudes. The relevant amplitudes are 

(j)^^{u) = (l)^{u) = 6u{l-u). (4.3) 

Now momentum configurations of n 7^ 1/2 contribute and the M^q*^" helicity amplitude is 
no longer zero. 

Figure ^ shows predictions for the cross-section. The M^q"' amplitude has a less steep 
t dependence than M^. The higher suppression of the chiral odd term is what we expect 



since in the definition of the lap integral in r space (see (2.7)), the chiral even amplitudes 
have a Ki function and the chiral odd ones a Kq. The Tqq element is the ratio of the 
longitudinal component of the cross-section to the total. Figure ^ demonstrates that this 



''Note that in the authors impUcitly equated the tensor vector meson couphng, fy, to the vector 
couphng, fv, rather than using the value estimated from QCD sum rules. 
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Figure 2: Collinear and 5-function approximation: da/dt for parameter values, {af^^^, m) = 
(0.2, mp/2). 
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Figure 3: Collinear and (5-function approximation: r-matrix elements. The predictions are param- 
eter independent. 

longitudinal component is numerically sub dominant, but becomes increasingly important 
with rising \t\. The change in the t dependence of the combined cross-section, due to this 
component, results in the previous fit of 7 = and fixed a^^ now being marginal. We now 
fit the data better by either running with t or a^^ with t. 

Note that the amplitude is still zero in our collinear scheme. This means that 
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Figure 4: Collinear approximation: da/dt for parameter values, {af^^^, m) = (0.2, mpjl). 

the r^_\ r-matrix element is still predicted to be zero. As observed for the case of the 5- 
function distribution amplitude, this is precisely the r-matrix most incompatible with such 
a prediction. Figure ^ shows the effect on the r-matrix elements of running the scale 
with i versus keeping it fixed. We observe that both provide a good fit for rgg and poor 
predictions for r^Q and r^^. The poor quality of fits to the r-matrix elements motivates 
the progression to higher twist. Higher twist is the only solution to producing a non-zero 
contribution. 

4.3 Higher twist effects 

Higher twist effects are introduced by relaxing the assumption that a collinear quark and 
antiquark exit the hard subprocess (i.e. we no longer have dipoles of zero transverse size) 
or by inclusion of higher Fock states (e.g. qqg) in the wave function. 

Potentially, we must consider the full six helicity amplitudes. The complete set of 
relevant asymptotic distribution amplitudes are 

4>\\{u) = 4>±{u) = 6uu 

g'f\u) =6nn ) (^) = ^(1 + ^ 

h{u) = l{l-e), (4.4) 

^Note that a flip in sign of either M+o or both M++ and M+_ would flip the sign of our prediction for 
rgo- We note that our two-gluon leading twist amphtudes agree with those of jl^] and we checked that our 
BFKL result reduces to the two-gluon prediction in the limit z 0. In addition we have checked that all 
our conventions agree with those of ZEUS, i.e. the direction of q, and that we have the correct phases for 
the polarisation vectors. 
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Figure 5: Collinear approximation: r-matrix elements for parameter values, {af^^^, m) 
(0.2, mp/2). 



where ^ = 2m — 1. Substituting these into (|2.lD "(pi^) we obtain the ampHtudes which were 
presented exphcitly in Q. We do not hst them again here. 

We want to explore the effect of higher twist terms on our observables by systemati- 
cally adding suppressed terms. Any observable is proportional to products of the helicity 
amplitudes i.e. MijAI^^. Consider the (parton level) cross-section for the production of a 
longitudinal meson: 



dcr+o 
dt 

(twist-2 X twist-2) (twist-2 x twist-3) (twist-2 x twist-3) (twist-3 x twist-3) 

In the expression above we see that the first term is leading (twist-2 x twist-2), the second 
and third l/|q| suppressed (twist-2 x twist-3), and the fourth l/|qp suppressed (twist-3 
X twist-3). Note that a (twist-2 x twist-4) term would also be l/|qp suppressed. Next- 
to-leading power behaviour in observables is therefore the highest level in which we are 
able to compute. Note that the helicity double-flip amplitude is subleading at this level of 
approximation as far as the total cross-section is concerned but that it provides the leading 
behaviour for r^^^. 

Figure ^ shows da /dt for the varying levels of approximation. The effect on the com- 
bined cross-section is predominantly that of normalisation. We see that adding in the 
(twist-2 X twist-3) to the (twist-2 x twist-2) terms has a significant effect, but that adding 
those next-to-next-to- leading terms that we have knowledge of has much less effect. This 
reassuring result suggests that the total cross-section is insensitive to the twist-four terms 
that we have neglected. 
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Figure 6: Higher twist effects: da/dt for parameter values, (af^jaf^^^, A^, m) = 
(0.17, 0.2, ml, mpjl). 

Figure ^ shows that the observables Tqq and are also relatively insensitive to next- 
to-next-to leading terms. The r^*!^ proves more sensitive. Adding the (twist-3 x twist-3) 
terms has an effect of order 100%. It is not surprising that this observable is the most 
sensitive to higher twist, since it is zero at leading twist. It therefore appears that we 
cannot trust our predictions for this r-matrix element and that it is not a good observable 
with which to test the validity of BFKL dynamics. However, Figure ^, makes it clearer 
what is happening. It shows the ratios of even to odd helicity amplitudes at rc = 0.1^. Our 
twist counting scheme ignores the effect of |q|-dependent effects in the hard scatter. The 
plots in Figure ^ are a test of the validity of this assumption. If effects of the hard scatter 
were the same, we would expect twist counting to give us M^^'^ /M'^q and M^^^ /M'^^ 
approximately oc \/— t and M'^_^'^ / M'^_^ oc l/\/— i. Our assumption seems reasonable for 
the (-1-0) and (-1 — ) amplitudes, however, for (-I--I-) the odd and even terms have a similar 
t-dependence. M^^^ is twist-3 but the q dependence from the hard scatter boosts its 
importance to the same significance as the twist-2 pieces. The thick line in the plot of 
Ti^i in Figure |^ is the prediction where we, in addition to the M^'^M^'^" term, kept the 
j^even j^fwen j^ov t significantly above IGeV^, we can see that the effect of terms 

genuinely l/|qp suppressed and above, is minimal. Thus, though rj*!^ is the most sensitive 
observable to higher twist effects, we can claim to have control over it. We shall therefore, 
from here on, use the formulae (|2.1| )- (p.6D with the distribution amplitudes computed up 

to twist-4 for M_| and up to twist-3 for all others. Our control over higher twist suggests 

that higher twist effects are unlikely to provide the solution to the fact we cannot fit the 
complete set of r-matrix elements. 

^At HERA, the mean x values that contribute for each amplitude are (x) ~ 0.1 



- 12 - 



0.2 

0.1 



-0.1 
0.1 



1-0.1 




-0.1 
'-0.2 
-0.3 





up to (twist-2 X twist-2) 
up to (twist-2 X twist-3) 






up to (twist-3 X twist-3) 

up to (twist-3 X twist-4) 


1 , 1 , 










1 


1,1," 




2 3 


4 5 ( 












1 


1,1, 




2 3 


4 5 ( 


1 1 1 1 1 1 1 1 1 






— up to (twist-2 X twist-3) + "extra" piece 








1 1 1 , 1 1 1 , 


1 


2 3 


4 5 f 



Itl [GeV] 



Figure 7: Higher twist effects: r-matrix elements for parameter values, (af^^^, A^, m) 
(0.2, ml, mp/2). 




Figure 8: Higher twist effects: even to odd amplitude ratios for parameter values, 
(af A2, m) = (0.2, m^., mp/2). 



4.4 The sensitivity of predictions to parameters 

We demonstrated in the previous section that, away from low i, our predictions are only 
weakly sensitive to higher twist corrections. The fit for the Tj'q r-matrix element remained 
particularly poor, however. In this section we explore the sensitivity of our results to vary- 
ing the various parameters at our disposal. We shall use a fixed a^J^ = 0.17, and it should 
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then be understood that only the shape of the i-distribution matters (its normaUsation 
being adjustable by varying a^^). Of course a^^ cancels in the r?^. 

4.4.1 Sensitivity to af^^^ 

Figure ^ shows a central curve corresponding to: {a^^ , af^^^, A^, m) = 

(0.17, 0.2, TTip, mp/2). The two other curves show the effect of varying the value of af^^^ 

by 25%. The effect is primarily that of normalisation. 



The r-matrix element predictions are shown in Figure IC. Tqq is the most sensitive 



to variations in af^^^, followed in order by r^Q and rj*!^. We find rjl^ to be relatively 
insensitive to this parameter whereas we remind the reader that it proved most sensitive 
to higher twist. We note that while the general quality of the fits for all the r-matrix 
elements is unaffected by varying af by 25%, a large value of Ug is associated with a 
larger suppression of the longitudinal cross-section. We also note that the fit to r^Q remains 
poor. 



4.4.2 Sensitivity to 7 

The default set, (al^ , af^^^, A^, m) = (0.17, 0.2, ml, mp/2), that provides the solid 



curve in Figure 11, is only marginally compatible with the cross-section data. This is 
because, as previously noted, the longitudinal contribution becomes increasingly important 
at higher t and this flattens the t dependence. The dotted curve on the same graph 
demonstrates that we can fit the data much better by running the BFKL scale, A 



2 



The predictions in Figure 12 for the r-matrix elements again show the sign problem 
in the prediction for r^^. Running A^ has not solved this. Note also that the predictions 
for rQQ and r^Q are sensitive to the values of the BFKL parameters and in particular, the 
quality of the fits becomes worse when A^ is run with t. We see that when A^ is run, the 
longitudinal fraction grows quicker than is compatible with data. 

4.4.3 Sensitivity to z rapidity 

Note that the rapidity variable z increases (or decreases) as we raise af^^^ (or 7). It is 
really this z variable that we have been investigating in the previous subsections. Figure 
[l3| contains a plot for each helicity amplitude at the rapidity values, z = 0.5, 0.75 and 1.0. 
The number in the brackets in the legends correspond to the areas underneath the curves. 
Note that the position of the peaks do not shift as we change rapidity. As z decreases 
all amplitudes fall, but a look at the integrated values demonstrates that the longitudinal 
amplitudes fall more slowly. If we lower z we raise the ratio of the longitudinal to the 
transverse contributions. 

This observation allows us to play off the effects of varying af^^^ and A^ simultane- 
ously. A non-zero 7 value in A^ improves the t dependence of the differential cross-section 
and in conjunction, increasing af^^^ drives the longitudinal component down. The fit 
values {al^ , af^^^, A^) = (0.17, 0.25, — t) therefore lead to an improved overall fit. 
The predictions for these values are shown in Figures 14-^. Let us here recall the familiar 



result that the dependence of all amplitudes on the partonic collision energy s is charac- 
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Figure 9: The sensitivity to varying af^^'": da/dt for parameter values, (af^, A^, to) = 
(0.17, ml-t, mp/2). 




Figure 10: Tire sensitivity to varying af^^^: r-matrix elements for parameter values, (A^, m) 
{ml - t, TOp/2). 



terised by the BFKL exponent, i.e. Mij ~ with A = 12 In 2 a^^^^/vr, provided that s 
is sufficiently large, say z > 1. 

4.4.4 Sensitivity to quark mass 

The sensitivity of our BFKL predictions to the quark mass parameter has not so far been 
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Figure 11: The sensitivity to varying 7: da/dt for parameter values, (af^, af^^^, m) 
(0.17, 0.2, mp/2). 
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Figure 12: The sensitivity to varying 7: r-matrix elements for parameter values, {af^^^, m) 
(0.2, mp/2). 



explored. We do so now. Figure 16 shows differential cross-section predictions. The central 
curve is that of the default set, where (a^'^, af^^^, A^, m) = (0.17, 0.2, m^, ■mp/2). The 
other two curves were generated by doubling and halving the constituent quark mass. 
Bearing in mind the fact that we have varied the mass parameter considerably, we see that 
the cross-section has a fairly robust prediction, its shape being fairly insensitive to changes 



-16- 






-0.25 
■ -0.5 
-0.75 
-1 















- z^l.O (-0.27) 






- (-0.11) 






.. 2=0.5 (-0.06) 




" , 1 , 1 , 


,1 





1 ' 








- z=1.0 (^0.29) 






- z=0.75 (-0.12) 






z=0.5 (-0.071) 






1,1. 










Figure 13: The sensitivity to varying z: the six hehcity amplitudes differential in u. The y-axis is 



dM 



du 



where Cp = isai'' {N^ - l/N^)eQpfp and (m, |q|) = {mp/2, VlOGeV) 




Figure 14: An improved fit: dci/di for parameter values, {a 
t, mp/2) 



IF BFKL 



A2, m) = (0.17, 0.25, m 



in m. 



Figure |17| shows the r-matrix element predictions for various quark masses. They are all 
sensitive to this parameter, both in shape and in normalisation. Figure |l^ helps to explain 
why this is so. It shows the six helicity amplitudes, differential in u, for m = rUp/A, mp/2 
and rup. We point out that it is not in fact the value of m that is significant, it is the 
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Figure 15: An improved fit: r-matrix elements for parameter values, (a^ , A , to) 
(0.25, ml-t, mp/2). 



value of m/|q|, and in these plots we fixed |q| = ylO GeV. Notice that as we increase 
the scale m/|q|, the peaks of the distributions become concentrated increasingly toward 
u = 1/2, as expected. However, the effect is quite dramatic. Halving the constituent 
mass strongly enhances the contributions close to the end-points. The effect is particularly 
marked for the longitudinal helicity amplitudes, since they are forced to zero at exactly 
u = 1/2. Raising the quark mass is therefore a mechanism for suppressing the longitudinal 
components with respect to the transverse. The dramatic dependence of the end-point 
region on the quark mass, especially for the M^^", and M^*^ amplitudes, is closely 

related to the observations of who found that these three amplitudes are divergent in 
the two-gluon exchange approximation due to divergent end-point behaviour. We here 
note that although our amplitudes are no longer divergent, even in the massless quark 
limit, they do retain a memory of these large end-point contributions. We explore the 
end-point region in much more detail in Appendix B. In passing we note that this physics 
is reminiscent of that for the DIS process 76 77*6 eff, where / is a fermion. The 
structure function is zero at u = 1/2 and has large contributions close to the end-points 
that become suppressed as we increase the fermion mass []T7| . 

It is interesting to note that for a large enough m, the M+o amplitude not only is small 
but also changes sign. Then, the prediction of r^^ flips sign to that of the data. Figure 



17 shows this begins to occur for m as low as rup, in the low |t| region. The evidence of 
this can also be seen even at large t, in Figure 18, where the values of M+o differential 



in u, become positive for small u. So, increasing m would be a possible mechanism for 
improving our fits, however, we find it difficult to justify a mass parameter so far away 
from the constituent mass. 
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Figure 16: The sensitivity to varying m: da/dt for parameter values, (af^, af^^^,A^) = 
(0.17, 0.2, ml). 

Increasing the mass has the effect of suppressing large dipole configurations due to the 
McDonald Bessel function, i^j(m|r|), present in the r space representation of the photon 
wavefunction (e.g. see (p.7|)) which has a large argument suppression ~ e"'"''"!. Figure 



18 illustrates the point; larger masses suppress the end-point contributions more strongly. 
End-point contributions in u correspond to large dipole sizes and the data may therefore 
be interpreted as suggesting that large dipole sizes should be more heavily suppressed 
than our calculations predict. The 5-function distribution amplitude that we used in the 
collinear approximation is an extreme way of achieving this suppression (by forcing u = 
1/2). We suggest that the real physics of this suppression arises as a result of Sudakov 
factors associated with the quark lines pgjj . 

4.5 Sensitivity to distribution amplitudes 

In the next part of our analysis, we examine the effect of subasymptotic corrections to the 



meson distribution amplitudes. We use the full distribution amplitudes of |12, 13 1 presented 
in Appendix ^ in conjunction with the helicity amplitudes ( |2.1| )-( |2^ . We remind the 
reader that we do not adjust any of the parameters of the distribution amplitudes. 

Figure ^ shows the differential cross-section in two distribution amplitude prescrip- 
tions which we label "asymptotic" and "BBKT" (Ball, Braun, Koike, Tanaka). The more 



accurate distribution amplitudes prove to make no qualitative difference. Figure 20, for 
the r matrix elements, show that the effects are quite modest in these ratios. 



Figure |2l| shows the six helicity amplitudes, differential in u, at fixed rapidity z = 0.75. 



Referring back to the amplitudes M^^^ and M^'"_f", given by (2^) and ( |2.3D , we see that 
the difference between their u dependence is a relative switch in sign between two pieces 
(one corresponding to the vector component of the Fierz decomposition and the other to 
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Figure 18: The sensitivity to varying to: the six hehcity amphtudes differential in u. The y-axis 
is ^^^^%r— ' where Gp = isal^'' [N"^ - l/N^)eQpfp and (z, |q|) = (0.75, VTOGeV). 



the pseudovector) . The plots in Figure ^ demonstrate a cancellation in the effects of the 
non- asymptotic distribution amplitudes for M^^", but an amplification for M^^"'. The 
effects of the different distribution amplitudes on the u distributions of the amplitudes are 
most significant for M^^'^ and M'^^. However, even then there is a significant cancellation 
in their integrated values. We note that subasymptotic corrections are no aid to improving 
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Figure 19: The effect sub-asymptotic distribution amplitudes da/dt for parameter values, 
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Figure 20: The effect of sub-asymptotic distribution amplitudes ||12|, 13 
parameter values, {af^^^, A^, m) = (0.25, — t, mp/2) 



-matrix elements for 



our fit to Tj'q. 

Before proceeding, we ought to comment that in our approach it is assumed that 
the dipoles which scatter elastically at high t are small in comparison to the meson size. 
Therefore the dependence of the vector meson wavefunction on the dipole size is neglected. 
On the other hand, we use the perturbative photon wavefunction, with a constituent quark 
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Figure 21: The effect of sub-asymptotic distribution amplitudes u% : the six heHcity amphtudes 



differential in u. The y-axis is ^ 
(0.75, mp/2, VWGeV). 



dM. 



du 



, where Cp — isa'/ {N^ — l/N^)eQpfp and (z, to, |q|) 



mass giving the upper cut-off of the dipole sizes in the photon (recall, that Ka{mr) — > 
exp(— mr) for mr ^ 1). Sensitivity of the amplitudes to this parameter turns out to be 
significant. This suggests that an additional suppression of larger dipoles by the transverse 
part of the vector meson wave function could well have a sizeable effect. 



5. Phenomenology for the p, and J/ip 

We now turn our attention to the (p &iid J/tp, reverting to the use of the asymptotic 
distribution amplitudes since wish to treat the three mesons on the same footing. Given 
the smallness of the sub-asymptotic corrections for the p (as illustrated in the previous 
section) this seems quite reasonable. 

To what extent do the observations, made for the p in the preceding sections, hold for 
the other mesons? The calculation for a general meson, V, with asymptotic distribution 
amplitude, requires only the substitution {fp/ fj, Qp, mp) — > {fv/ fv^ Qv, my). The val- 
ues for the meson decay constant (fv) and electromagnetic coupling {Qv) are constants of 
proportionality, my appears in the BFKL logarithm and the 'jV impact factor (through 
the quark mass m = mv/2), and is qualitatively the most significant parameter. The mass 
of the (p is similar to that of the p, while the mass of the J/V' is significantly larger. We 
might therefore expect the predictions for the (p to be, qualitatively, similar to the p, and 
those for the J/ip to be, perhaps, very different. This is in fact what we shall see. 



Figures 22 and 23 reproduce our 'improved fit' prediction for the p shown in Figure 
14, where we played off the effects or adjusting af^^^ and to our advantage. Now 
we also show the two-gluon exchange predictions. Two-gluon exchange with a fixed strong 
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coupling predicts a differential cross-section far too flat in Running the coupling solves 
this problem and provides a good fit. However, looking at the two-gluon predictions for 
the r-matrix elements (which are independent of how we treat ctl^), we see that the values 
for the Tqq far exceed those constrained by the data. This means that the longitudinal 
component dominates for two-gluon exchange. This is in line with our observation that the 
longitudinal fraction increases as we lower the z rapidity (recall that the BFKL solution 



tends to two-gluon exchange in the limit z ^ 0) and is illustrated explicitly in Figure 24 



Thus, even though we have succeeded in getting the two-gluon curve to agree with the 
data for da/dt it is generated by fundamentally wrong dynamics. We do however note 
the improvement with respect to working within the collinear and 5-function distribution 
amplitude approximation. In that prescription, two-gluon exchange was studied in |, ||. 
The dip to zero manifest at \t\ = my is not seen in data and is a somehwat artificial 
prediction arising from the crudity of the treatment of the meson distribution amplitudes. 



Figures |25H26| and |271-Pq present our results for the (f) and J mesons. The BFKL 
curves which we label "(1)" in Figures ^ and ^ were obtained by simply making the 
appropriate change in meson constants and masses. No further adjustment of the other 
parameters is made relative to the fit for the p meson. For the curves labelled "BFKL (2)" 
we altered the values of the meson tensor coupling but have kept the ratio 



= 'i-fv. (5.1) 

Jp 

Fixing fy in this way improves the fit to data. Note that if we allowed ourselves slightly 
more freedom with the values of fy, the quality of the fits to the t-distribution could be 
improved still further. 

The (p predictions are similar to those for the p. We again see that the two-gluon 
exchange t-distribution requires a running strong coupling to fit the data and that the r- 
matrix elements point to this being an accident. The BFKL predictions are again correctly 
dominated by the transverse contributions, but again fail to predict the correct sign for 
r^Q in the case of the (j). However, for the J/V'j the BFKL predictions are compatible 
with all the observables. The large quark mass drives the longitudinal amplitudes small 
enough to agree with the Tj'q data. We note that although the two-gluon exchange r-matrix 
predictions for the J/ijj are in marginal agreement with the data, we have to doubt the 
validity of the underlying dynamics, given the comments above. 

Comparing two-gluon exchange and BFKL predictions, we see that the BFKL calcu- 
lation is a definite step forward. The success of the our predictions for the (j) and J /ijj \s a 
noteworthy result. We now understand the success of the apparently naive analysis of Q 
as being due to the fact that it correctly identified the dominance of the M'^^ amplitude. 
However, the inability of BFKL to agree with the data for Tj'q suggests that while the 
longitudinal contribution is brought down by BFKL effects, it is still not under enough 
control to accurately describe all the observables measured at ZEUS and HI. 
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Figure 22: p photoproduction: da/dt. The fixed two-gluon curve was calculated 
with the parameter values (af^, /p, /J, m) = (0.27, 0.216 GeV, 0.160 GeV, mp/2). The 
run two-gluon curve was calculated with the parameter values (af^(l GeV), /p, /J, m) 
=(0.30, 0.216 GeV, 0.160 GeV, TOp/2). The BFKL calculated for the parameter values 
{ai^, af^^^, A^, /p, /J, to) = (0.17, 0.25, - 0.216 GeV, 0.160 GeV, TOp/2). 



6. Summary 

We have compared our theoretical predictions for da/dt and the three r-matrix elements 
to the data of ZEUS and HI, in various levels of approximation. The only meson we found 
complete agreement for, in any scheme, was the J/tp. We nevertheless could generally 
obtain good fits to the Tqq observable. Fits to the cross-sections favoured running the 
BFKL scale in a natural way: = - t. The theoretical predictions for r^^ were 
the main obstacles to obtaining satisfactory fits for the p and (p. We found the natural 
predictions of both two-gluon and BFKL pomeron exchange to be the wrong sign. The 
predictions of the observables proved insensitive to higher twist terms, suggesting that these 
corrections are unlikely to provide a solution. We also found that using the distribution 



amplitudes of |12|, |13[] had little effect on the quality of fits compared to the asymptotic 
ones. The inability of theoretical predictions to fit r^Q, and the general sensitivity of the 
r-matrix elements to the quark mass parameter, m, therefore seem to locate the problem 
at the level of the hard subprocess. We suggest that the data hint that we may require 
a greater suppression of larger dipoles than we predict and that this may be provided by 
Sudakov suppression of emissions off the quark lines. 

We note the dominance of the M°'^ contribution and have demonstrated that this 
corresponds to the amplitude used in the BFKL calculation of |^ (which was, in addition, 
conducted in a (5- function approximation). We were able to get excellent agreement for 
the 4> and J/ip differential cross-sections using the same parameters that fitted the p. On 
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Figure 23: p photoproduction: r-matrix elements. The two-gluon curve was calculated with the 
parameter values (/p, /J, m) = (0.216 GeV, 0.160 GeV, mp/2). The BFKL curve was calculated 
for the parameter values (af^^^, A^, /p, /J, m) = (0.25, - t, 0.216 GeV, 0.160 GeV, TOp/2). 
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Figure 24: The relative contributions for two-gluon exchange and full BFKL (z = 1.2, 

t = —10 GeV^): the six helicity amplitudes differential in u. Note that the two-gluon exchange 

results have been multiplied by a factor of 3. 



the whole the evidence suggests that leading logarithm BFKL predictions work well for 
the transverse amplitudes but, despite suppressing the longitudinal component relative to 
Born level, it fails to have enough control over this contribution. 
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Figure 25: 4> photoproduction: da/dt for parameter values. The fixed two-gluon curve was 
calculated with the parameter values (af^, Z^, /J, m) = (0.25, 0.231 GeV, 0.215 GeV, 
The run two-gluon curve was calculated with the parameter values (^^^(IGeV), Z^, J^", m) 
= (0.28, 0.231 GeV, 0.215 GeV, TO0/2). BFKL prediction (1) is for the parameter values 
(af^, af^^^, A^, /J, to) = (0.17, 0.25, m^-t, 0.231 GeV, 0.215 GeV, TO0/2). BFKL predic- 
tion (2) is for the parameter values (af-^, af A^, /J, to) = (0.17, 0.25, -f, 0.231 GeV, 
0.171 GeV, TO0/2). 
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Table 4: Classification of distribution amplitudes up to twist-4. 
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A. Vector meson distribution amplitudes 



The authors of presented results for the complete set of distribution amplitudes up to 
twist-3 for the p and (f> (and K) [12|. These results were extended to twist-4 in [^]. The 
relevant distribution amplitudes, up to and including twist-4, are classified in Table ^ 




Figure 26: cj) photoproduction: r-matrix element. The fixed two-gluon curve was calculated with 
the parameter values (J^, fj, m) = (0.231 GeV, 0.215 GeV, m0/2). The BFKL prediction ((2)) is 
for the parameter values (af A^, /J, m) = (0.25, ml - t, 0.231 GeV, 0.171 GeV, m^/2). 





a}\ 


«| 


Cs 




^3. 


p 





0.18 ±0.10 


0.032 ±0.010 


-2.1 ± 1.0 


3.8 ± 1.8 



Table 5: Chiral even parameters: /i^ = 1 GeV^. 







ai 


Cs 


T 






p 





0.20 ±0.10 


0.032 ±0.010 


7.0 ±7.0 


0.10 ±0.05 


-0.10 ±0.05 



Table 6: Chiral odd parameters: fi^ ^ 1 GeV^. 



We quote all parameters at the renormalisation scale fi^ = 1 GeV^ and neglect the 
slow running. In the q ^ oo limit, the parameters vanish. We refer to this as the asymp- 
totic limit, and the distribution amplitudes in this limit as the asymptotic distribution 
amplitudes. 

Tables |5| and |6| hold the values of all the relevant parameters, at the scale /i^ = 1 GeV^. 
We do not include quark mass corrections in this study. They are in any case zero for the 
p meson. 
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Figure 27: J/-0 photoproduction: da/dt. The fixed two-gluon curve was calculated with the pa- 
rameter values (af^, /j/^, /J/^, w) = (0.23, 0.405 GeV, 0.405 GeV, mj/^/2). The run two-gluon 
curve was calculated with the parameter values (^^^(IGeV) /j/^, /J/^, m) = (0.29, 0.405 GeV, 
0.405 GeV, m.7/^,/2). BFKL prediction (1) is for the parameter values {a[^,a^^^^, A^, 
/j/^, /J/^, f") = (0.17, 0.25, m?ji^ - t, 0.405 GeV, 0.405 GeV, mj/^/2). BFKL prediction (2) is 
for the parameter values (af^, af"^^, A^, /j/^, /J^^, m) = (0.17, 0.25, m]^^ - t, 0.405 GeV, 
0.300 GeV, TOj/^/2). 

A.l Twist-2 distribution amplitudes 
A. 1.1 Chiral Even 

There is one chiral even, twist-2 amplitude: 

<A||(n) = Quu {l + Safe + a| ^(5^^ _ _ (^.1) 
and the asymptotic distribution is 

(/)J'^(u) = Quu. (A.2) 

A.l. 2 Chiral Odd 

There is one chiral odd, twist-2 amplitude. This can be written 

(/)_L(n) = Quu + Safe + ^(5^^ _ _ (^.3) 
The asymptotic distribution is 

(j,'''y{u) = Quu. (A.4) 
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Figure 28: J/i/j photoproduction: r-matrix elements. The two-gluon curve was calculated with 
the parameter values {fj/-4,,fj/^,m) = (0.405 GeV, 0.405 GeV, toj/^/2). The BFKL predic- 
tion ((2)) is for the parameter values {af^'^'^, A^, fj/^, fj^^, m) = (0.25, m^^^ - t, 0.405 GeV, 
0.300 GeV, mj/^/2). 



A. 2 Twist-3 distribution amplitudes 
A. 2.1 Chiral Even 

There are two even, twist-3 distributions: 



g^^\u) = Quu 



1 + <! -4 + -(3(1- — + — ^•^ 
'42^3^1 16 ^ 16 ^ 



{5e - 1) 



(A.5) 



and 



5?(^) = |(l + a + (fa| + 5C3) (3^2-1) 
15 



+ 



112 



"2 + ^^3(3w|' 



(3 - 30^^ + 35^ 



(A.6) 



where the subscript of '3' refers to three-particle corrections and we neglect mass correc- 
tions. In the asymptotic limit 



(A.7) 



and 



(A. 
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A.2.2 Chiral Odd 

There is one relevant chiral odd twist-3 distribution: 

h\;^ (u) = 3e + l fiBf - 3) + ^ Cs {3 - 30e + 35^). (A.9) 
In the asymptotic limit 

h'^^'''y{u)=3f. (A. 10) 

A. 3 Twist-4 distribution amplitudes 
A. 3.1 Chiral Odd 

There is one relevant chiral odd twist-4 distribution; 



h^iu) = 1+^-1 + -at - mi + Ci ) j C.'^iC) 

+ (-^4-^C3-j)cr(0. (A.11) 
where C^(^) are Gegenbauer polynomials. In the asymptotic limit 

hi^y = i-cl/\o = l{i-e)- (A.12) 

A. 4 Distribution amplitude schemes 

We have six helicity amplitudes to consider, i.e. ( |2.1| )-( ^l6| ). Their individual dependence 
on the distribution amplitudes is as follows: 

nTW^ (^^^-(l-2n)j[|"dt;(<^||W-<7?W)^ (A.14) 

^e.en(^)^ (^^^ + (1 " 2n) ^" dt; (</)|| (t;) - 5? (t;))^ (A15) 

$^o^(n) ^ £dv (^h'^\v) - ^^{v)) (A.16) 

n'ti^) ^ Mn) (A.17) 

^f'(n) ^ ^dvj^ri (^h\f\v) - 1<P±{V) - ^/i3(r/)) • (A.18) 
We now present the explicit formulae for these amplitudes in four different prescriptions. 
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A. 4.1 Prescription 1: leading twist (collinear) approximation with (5- function 
distributions 

We keep only the twist-2 contributions and equate the corresponding distribution amph- 
tudes to a (^-function that enforces the quark and antiquark share the meson momentum 
equally. Then 



^eyen(„) = (A.21) 
\odd 



^"^{u) = (A.22) 



K'ti^) = 6(u-\] (A.23) 



^°^{u) = 0. (A.24) 

A. 4.2 Prescription 2: leading twist (collinear) approximation with asymptotic 
distribution amplitudes 

We keep only the twist-2 contributions and equate the corresponding distribution ampli- 
tudes to their asymptotic forms. Then 

= Quu (A.25) 

^l"!!\u) = (A.27) 

$fo'^(tx) = (A.28) 

= 6uu (A.29) 

(«) = 0. (A.30) 
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A. 4. 3 Prescription 3: higher twist approximation with asymptotic distribution 
amplitudes 

We keep all twist and equate the corresponding distribution amplitudes to the asymptotic 
forms. After integration: 

= 6nn (A.31) 
= 3uu{u^ + n^) (A.32) 
= 6u^u^ (A.33) 
$f(f (u) = 3(1 - 2u)uu (A.34) 

(-u) = ^u\^ (A.35) 

A. 4. 4 Prescription 4: higher twist approximation with BBKT distribution am- 

plitudes 

We keep all twist and equate the corresponding distribution amplitudes to the BBKT forms 
listed in A1-A3. We neglect to explicitly quote the trivial, but bulky, integrated forms. 

B. End-point behaviour 

B.l Two-gluon exchange and Reggeisation corrections 

The importance of contributions at, and close to, end-points is a matter of debate (see for 
example (|6|, ^, ^])). The authors of ^] worked in the two-gluon exchange approximation in 
the light-quark hmit^^ and concluded that M^Y^", Mfff and M^*^ were all plagued by end- 
point divergences. They argued that these end-point divergences would be brought under 
control by Sudakov corrections and hence that the dominant amplitude would remain M^_^. 
In 0, the end-point divergences were interpreted as the signal that QCD factorisation was 



breaking down. In the context of a specific model of the vector meson wavefunction, |23| 
obtain a meson wavefunction, dependent on both longitudinal and transverse coordinates, 
by a Borel transform of the photon wavefunction. The model of [^3| gives an explicit form 
for the end-point contributions at all twist, it is shown that those contributions can be 
resummed and that the breakdown of factorisation by end-point divergencies at leading 
twist is only apparent. Here we discuss the role of end-points within the context of leading 
logarithmic BFKL resummation. 

In Figure we also observe the appearance of large end-point contributions (in the 
same amplitudes noted by Q). It is the case that these contributions lead to divergent 
matrix elements in the two-gluon approximation in the case of massless quarks. However, 



^We refer to our previous paper pi] for a discussion on this limit. 
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Figure 29: The M™^" amplitude against to, the quark mass parameter. The fixed parameter 
values are {z,mv, \t\) — (1.0, mp, lOGeV^). The constant, Cq = sal{N^ — 1/N^)fvmveeq and the 
asymptotic distribution amplitudes were used. 

these divergences are not present after re-summing the BFKL logarithms. In Figure 29 
we show that the the BFKL prediction for the (++) amplitude remains finite even in the 
massless quark limit. We shall investigate how this divergence in (++) at the Born level 
arises and how higher orders bring it under control. The authors of worked with the 
asymptotic distribution amplitudes; we shall also do so (see Appendix |A| for the asymptotic 
formulae). The end-point behaviour is independent of this approximation. 

The light-quark limit corresponds to putting the mass of the quarks that form the 
meson to zero. This leads to the replacements 

mKi{m\r\) and Ko{m\r\) —mln{m\r\) ~ 0, (B.l) 

where we have assumed that m\r\ — > 0, i.e. that we get no significant contributions from 
asymptotically large transverse sizes. 

We can apply the light quark approximation to the Born (two-gluon) level amplitudes 
previously derived. For M'^^"''^^ one has 

^e.en2, ^ _^CvfvMv^, (in (l^J^) _ 1 + 2n„ J . (B.2) 

A clear divergence as the cut-off Umin 0. 

The full BFKL solution has proven in the past to improve infra-red finiteness of scat- 
tering amplitudes; see for example ||2^, ^ for the case of qq qq scattering. At any 
finite order in this scattering amplitude diverges. Even in infra-red finite observables 
the Born level result can be seen to give undue weight to contributions from momentum 
configurations which are unphysical, i.e. where all the momentum is short-circuited down 
one gluon. 
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The authors of |^] demonstrated how this problem can be resolved without recourse 
to the full BFKL result, but rather by a fixed order resummation of a subset of the BFKL 
logarithms. Following Q, we go from the Born amplitude to the resummed one by the 
substitution 

k2(k-q)2 ^ k2(k-q)2 V J ' 

The substitution above clearly suppresses small gluon momenta and the formerly divergent 
amplitude can now be written (using the definition of the incomplete beta function) 

j^^ven2gR ^ _^^CvfvMv 

1 Umin 

B{l-Umin,2z,l + 2z) - B{l-Umin,l + 2z,l + 2z) . (B.4) 

We can safely put Umin ^ for z > 0: 

^even29R ^ ^^^2z, I + 2z) - B{1 + 2z,l + 2z)) 

|qp 

evr T{2z)T{2 + 2z) 

= -J^.sCyfyMy ^ • (B.5) 

As it should be, we recover the divergence in the limit z O'J 

Ml-^^^9R^^^Q^^_^^CvfvMv-. (B.6) 

1^1 ^ 

We observe that this is equal to that derived for qq — > qq scattering amplitude in |2C, 21 1, 
excepting a factor of —3CvfvMv/2q'^. This is understood since in this limit the quark and 
antiquark are far apart and the dominant contribution comes from coupling the exchanged 
reggeised gluons to a single parton [^, The clear result is that the resummed (++) 
even amplitude remains finite even in the massless quark limit for non-zero rapidity. We 
have demonstrated analytically that higher orders accomplish this and so it is no longer 



surprising that Figure 29 demonstrates convergence at m — > 0. 
B.2 BFKL 

The previous subsection is sufficient to illustrate that the BFKL result is free from infrared 
divergences. Here we take a closer look at the end-point behaviour of the full BFKL 
amplitude. The reader not interested in the details may wish to jump directly to the 
final result, (|B.23| ). Equation (2.1) and the definition of the asymptotic distributions in 



Appendix |^ give us the result 

^e.en ^ ■!CvfvMv T rf^ 3 (^2 + -2) 







E 

i=— 00 



00 



, J" +n exp[x2riHf] r r o i\ 

[w^ + [n — 1/ lyWv^ + [n + 1/2) J smynrv) 



(B.7) 



^Of course we trivially reproduce the divergence of (|B.2|) if we first take the limit z — > 
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where 



r-C'+ioo J A 

lapiiy, n, q,u;a) = - :^r{a/2 - C)r(-a/2 - C) t< {i sign (1 - 2tz))«-^+" 

2 Jc'-ioo ^TTZ 



1^1 



C'-ioo 

4 



X { ) [sin7r(Q; + /x + C) B(a, n,q* ,u,C) B{P, n,q,u* ,C) 



-i-ir sin7r(a - + C) ^(a, "M, C) "M, Q, u\ Q ] (B.8) 



and 



with Tq = 4m^/|gp. 

The even amplitudes have a = 1 so the integrand in Iai3{i^,n,q,u;l) has a pole at 
1^ = —1/2, where — 1 < C < — 1/1. If we shift the contour to the right we pick up pole 
contributions from Re(^ > —1/2. In the limit m ^ we might expect the leading pole to be 
that at C = —1/2, the residue of which is independent of m, the other poles contributions 
vanishing in the limit. The integrand is complicated, however, expressed how it is in terms 
of blocks of hyper geometric functions. It is difficult to see clearly the analytic behaviour 
in For example the true expansion parameter may not be Am? /q^, rather something 
like Am? /v?q^, where it is not clear whether the 'leading pole' actually is leading in the 
simultaneous limit m — > and u ^ 0. We shall see if some manipulation of the formula 
can be enlightening. Using the following hypergeometric identity 

2^1 (a, 6; c; z) = li^ffi^ 2Fi{a,a + 1 - c;a + 1 - b; 1/z) 
i (oji (c — a) 

+ r(a)r(c-b) ^~^^"' 2F,{b,b + l-c;b + l-a; 1/z), (B.IO) 

we can re- write, 

B{a, fi, q*,u, = (-4uu)-('*+2+a+C)/2 j r(/x + 2 + a + C) 
r(-3/2-a-C) / 1 \ ^ 



Lr 



( 2 ' 2 , 5/2 + a + C, 4nu 

^ r(3/2 + Q + c) ( 1 V""''""^ 



4ltlt 



X2F1 ^ ^, ^,-l/2-a-C,4ti« 



(B.ll) 
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The advantage of this expression is that the hyper geometries now have well defined power 
series expansions for < u < 1. We can simplify this expression further using the gamma 
function identity, 



1 



22^- 



1 



r(z)r(z + i/2) r(2z)' 



(B.12) 



We find, 



TT 



B{a,i^,q*,u,C) ■- 

X [c(/x, a, |n|, + (-uu)-3/2-«-CL.(/i, a, \u\,0, 



(B.13) 



where we define 



C{n,a, \u\,C) 
X2F1 



r(/x + 2 + g + C)r(-3/2 -a-O 

r(/x-i-a-c) 

IJ, + 2 + a + C 2 + a + C-^ 



,5/2 + a + C,4:Uu] 



(B.14) 



and 



D{fi,a,\u\,C) =r(3/2 + a + C) 

>-1-q;-C l + a + C + n 

X2-ri 



, -1/2 - a - C,4'U'u 



(B.15) 



Note that both the C and D functions are single-valued in u, but that the B function 
explicitly breaks single- valuedness. This is not a problem, since it is only the combination 
of blocks that must be single- valued. The blocks occur in the combination 

sin7r(Q; + + C) B{a, fi, q*,uC)B{l3, /i, q, u* , () 

- (-I)"sin7r(a - + C) B{a, -/x, q*,uC)B{P, -fi, q, u*,C) 

C{ii,a, \u\,QC{fi,P, \u\,C) 

+C{„,a, \u\,ODiil,P, \u\,0 iuu)-'/W,±^m+PH) 
+Difi,a, \u\,OC{fi,P, \u\,0 (^x^)-3/2-"-CeT-(3/2+"+C) 

+D{f,, a, \u\, QDifl, /?, {uur^~--^-^H-^r~^ 



(-l)''sin7r(a!-/x-|-C) 



2\^1 



2-4-2C 



q J TT 



(B.16) 



This expression is complicated but in fact simplifies. Note that in the square bracket there 
are two multi-valued pieces corresponding to the cross multiplication of C's and D^s. The 
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sum of the two multi-valued pieces are not single-valued together in general; they have 
different n-dependences. In fact each multi-valued piece cancels exactly with another piece 
from the other square bracket®. We can throw away the cross multiplied terms. We obtain 
the result, 



Io,i3ii^,n,q,u;a) = ^ (isgn (1 - 2^))" 

pC'+ioo JA 

X / ^r(a/2-C)r(-a/2-C)(mVklY 
sin7r(a + /i + C) G{a, (3, fi, Jl, \u\, () 



-1)" sin7r(a — + C) G{a, (3, — /i, — /i, \u\, () 



(B.17) 



where we introduce the notation, 



G{a,f3,iJ,,fi, \u\X) = C{n,a, \u\,C)C{^, (3, \u\,C) 

+{-ir-f (nn)-3"°-/^-2C ^(^^ \u\,OD{ll, P, \u\,C). (B.18) 

We are now in a better position to examine the behaviour of the (++) even amplitude as 

Note that the full (++) even amplitude is proportional to the integral over u of /qo 
multiplied by the tt-dependence arising from the distribution amplitudes. At the end-points, 
we get the behaviour, 

±^ (X 3uu(u^ + u^) u-^~^^ 3n-2-2C_ g 

du 

The leading pole approximation would put C ~^ —0.5, leaving the integrated amplitude 
logarithmically divergent as in the Born level case. However, earlier we showed that re- 
summing leads to a finite result for non-zero rapidity. We now seek to show that the 
leading pole divergence, for non-zero rapidity, is an artifact indicating the break down of 
assumptions made rather than genuine asymptotic behaviour. 

The G-function (B.18) simplifies greatly at the end-points. In the limit u ^ the 



dominant u-dependence allows us to write 

G{a, f3, 12, Jl, \u\ ^ 0, C) ^ (-1)"-^ ^-3-a-/3-2C r(3/2 + a + C)r(3/2 + (3 + C), 

(B.20) 



^This has been verified numerically. 
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where we have put the hypergeometrics to unity in this hmit. This simphfication leads to 

2+/3 / r, \ 2+a 



C^^u, n, u;a) = ^ (^) (1) " (-1)^-^ (i sgn (1 - 2n))"-^- |^ 

f-C'+icO 

X / ^r(a/2-C)r(-a/2-C)r(3/2 + a + C)r(3/2 + /? + C) 



xsin7r(a + 1/2 + C) (mVkl^)^ 



where we have used the integral definition of the Meijer G-function For (a, a, (3) = (1, 0, 0) 
we find that this simplifies; 

Gil ( -^\'}'''-}''\ = G]] ( -^\7}A = . ^ .. . (B.22) 



'22 \ ^2|g|2llA-l/2 ) ^11 V^,2|^|2ll/2 



1 + 



We can therefore write down the (++) even amplitude in the u — > limit as 
dMf^'^{u 0) 3Cy/yMym2 ( 2^ I 1 

1 + 



n=+oo „oo 2 I 2 



X 

n=— oo 



y / dv—^ -— ^^^^ —:rre^^-^^>. (B.23) 

Note that the u-dependence has completely factorised from the n and v dependences. In 
fact the result, up to a factor, is that of qq qq scattering found in [20, |2l[|. We observed 
this earlier when performing the resummed calculation. From inspection of ( B.23 ), for any 
particular conformal spin, we can see that it vanishes as n — > for finite m/\q\ and also 
for m/\q\ and non-zero u. We know, however, that we must recover the Born level 
divergence as z — > 0. In fact we do achieve this; the sum over conformal spins is infinite 
for zero z. The fact that conformal spin has decoupled from u in our formula seems to 
imply that for zero z we face a divergence regardless of m ^ 0. In fact this is not so. For 
sufficiently large n the assumption we made in putting the hypergeometric functions to 
zero breaks down. We should have a cut off in n relating to the inverse of Umin- 

B.3 Conclusions 

In this appendix we have shown that 

• the two-gluon exchange (H — h) even amplitude diverges as Inumm as Umin — for 
m = 

• the resummed (++) even amplitude converges as Umin and m ^ for non-zero 
z, but diverges as l/z, as z — > (the two-gluon limit). 
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• the BFKL (++) even amplitude converges as Umin and m — > 0, when z is 
non-zero. 

• for the complete (or resummed) leading log calculation we can only get a divergence 
with the following simultaneous limits {umin, m, z) ^ and any one can provide a 
cut-off. 
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